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ABSTRACT 

We present an analysis of 18 short-cadence (SC) transit lightcurves of TrES-2b using quarter (QO) 
and quarter 1 (Ql) from the Kepler Mission. The photometry is of unprecedented precision, 237 ppm 
per minute, allowing for the most accurate determination of the transit parameters yet obtained for 
this system. Global fits of the transit photometry, radial velocities and known transit times are used 
to obtain a self-consistent set of refined parameters for this system, including updated stellar and 
planetary parameters. Special attention is paid to fitting for limb darkening and eccentricity. We 
place an upper limit on the occultation depth to be < 72.9 ppm to 3-cr confidence, indicating TrES-2b 
has the lowest determined geometric albedo for an exoplanet, of Ag < 0.146. 

We also produce a transit timing analysis using Kepler^s short-cadence data and demonstrate ex- 
ceptional timing precision at the level of a few seconds for each transit event. With 18 fully-sampled 
transits at such high precision, we are able to produce stringent constraints on the presence of per- 
turbing planets, Trojans and extrasolar moons. We introduce the novel use of control data to identify 
phasing effects. We also exclude the previously proposed hypotheses of short-period TTV and addi- 
tional transits but find the hypothesis of long-term inclination change is neither supported nor refuted 
by our analysis. 

Subject headings: planetary systems — stars: individual (TrES-2b) techniques: spectroscopic, photo- 
metric 



1. INTRODUCTION 

TrES-2b is a transiting planet discovered by the 
Trans-atlantic Exoplanet Survey (TrES) which happens 
to reside in the field-of-view for t he Kepler Mission^ 
(jBasri et all 120051 : iKoch et al.ll2007D. The fact that the 
planet was discovered by TrES (|0 'Donovan et al.ll2007t ) 
provides several advantages for Kepler. Firstly, the star 
was selected as one of the 512 targets for immediate 
short-cadence (SC) observations. Secondly, the planet's 
ephemeris is well-characterized from ground-based mea- 
surements meaning a search for long-term transit time 
variations (TTV) is possible. Thirdly, TrES targeted 
brighter stars than Kepler and thus TrES-2 is somewhat 
brighter (V=11.4) than typical Kepler stars (V = 12 to 
14). 

In iGilliland et al.l (|2010l ). a presentation of the first 
TrES-2b lightcurves was provided, but the focus of the 
paper was to demonstrate the properties of the SC data 
rather than a detailed study of the planet's properties. 
In this paper, we present a comprehensive analysis of the 
first 18 transits observed by Kepler in quarter (QO) and 
quarter 1 (Ql) in short-cadence mode. The photometry 
is analyzed in combination with the radial velocity (RV) 
data and known transit times of the system, and the com- 
bined results a llow for a refined YY-isochrone analysis 
()Yi et al.ll2001|) to derive a complete and self-consistent 
set of system parameters. Particular attention is paid 
to fitting for both eccentricity and limb darkening co- 
efficients, making the results as model-independent as 
possible. 



TrES-2 is somewhat remarkable, if nothing else, for 
having been the su bject of numerous tentative detec- 
tions. For example, iRaetz et al.l ()2009l ) claimed to have 
detected repeated dips in the lightcurve 1-2 hours af- 
ter the main transit event and pro posed a seco n d res - 
onant planet as an explanation. iRabus et al.l (|2009l ) 
claimed to have detected transit timing variations for 
TrES-2b of period 0.21 cycles (where one cycle is one or- 
bital period of TrES-2b) and 50 s amplitude and pro- 
posed a 52Mffl exom oon as a possible explanation. Fi- 
nally, iMislis fc Schmittl (j 20091) claim to have detected 
long-term inclination change in the system. In this work, 
we will also investigate the compatibility of these claims 
with the Kepler photometry. 

The SC mode was made available for the purposes 
of studying asteroseismo logy and transit ti ming varia- 
tions (TTV). Although iKipping fc BakosI (|2010l ) have 
shown that even the long-cadence (LC) is capable of per- 
forming TTV at the level of ~ 20 s, the SC data has 
the potential for further improvement on this. Such 
precision would allow for the detection of satellites 
(jSartoretti fc Schneideij|1999l:lKippineil2009D. Mars -mass 



pertu rbing planets ( Agol et al. 2005t Holman fc Murray! 
120051 ) and Trojan bodies (|Ford fc Holmanll2007l ). 
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2. DATA HANDLINC 

In this section, we will list the sequential steps we took 
in processing the Kepler photometry for TrES-2b. 

2.1. Data Acquisition 

We make use of the "Data Release 5" (DR5) from 
the Kepler Mission, which consists of quarter (QO) 
and quarter 1 (Ql). Full details on the data processing 
pipeline can be found in the DR5 handbook. Numer- 
ous improvements have been made over the previously 
available MAST (Multimission Archive at STScI) data. 
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including most relevant for this study an inclusion of B JD 
(Barycentric Julian Date) time stamps for each flux mea- 
surement. The previous version of the data only included 
cadence numbers and thus the inclusion of barycentric 
corrected time stamps is a marked improvement^. 

2.2. Long-Term Detrending with a Cosine Filter 

We make use of the "raw" (labelled as 
"AP_RAW_FLUX" in the header) data processed 
by the DR5 pipeline and a detailed description can 
be found in the accompanying release notes. The 
"raw" data has been processed using PA (Photometric 
Analysis), which includes cleaning of cosmic ray hits, 
Argabrightenings, removal of background flux, aperture 
photometry and computation of centroid positions. 

The data release also includes corrected fluxes (labelled 
as "AP_CORR_FLUX" in the header), which are out- 
putted from the PDC (Pre-search Data Conditioning) al- 
gorithm developed by the DAWG (Data Analysis Work- 
ing Group). As detailed in DR5, this data is not recom- 
mended for scientific use, owing to, in part, the potential 
for under/over- fitting of the systematic effects. 

For the sake of brevity, we do not reproduce the details 
of the PA and PDC steps here, but direct those interested 
to lGilliland et al.l (|2010f ) and the DR5 handbook. 

The QO and Ql PA photometry are shown in Figure [T] 
respectively. One challenge in attempting a correction 
is assessing which components are astrophysical in na- 
ture and which are instrumental. In general, we wish 
to preserve the astrophysical signal as much as possible. 
However, in practice, any signals occurring on timescales 
greater than the orbital period of the transiting planet, 
whether instrumental or astrophysical, have a negligi- 
ble impact on the morphology of the transit lightcurve, 
which is ultimately what we are interested in for this 
study. These signals may be removed by applying a high- 
pa ss filter to the photometry, in a similar way as was used 
bv lMazeh et al.l ()201Q) for the spaced-based CoRoT pho- 
tometry. 

To remove the long-term trend, vi sible in Figure 111 we 
applied a discrete cosine transform () Ahmed et al.lll974l) 
adopted to the unevenly spaced data. 

We first removed the 18 transit events with a margin 
of 6559.4s either side of the times of transit minimum. 
This value was chosen as it represents the lst-to-4th con- 
tact duration and thus we essentially remove the transit 
plus one half of the total transit duration either side of 
each event. We also remove outliers, identified as those 
points lying S-u away from a spline-interpolated running 
median of window-size 20 minutes. Treating QO and Ql 
separately, we fitted the remaining data with a linear 
combination of the first N low- frequency cosine functions: 

/.fe)=cos(^) (1) 

Where tj is the timing of the j"^ measurement, « = 0, 
in integer steps and N is equal to the rounded integer 
value of {2D/APp) where D is the timespan of the ob- 
servations and Pp is the orbital period of TrES-2b. For 
QO, we used Af = 2 and for Ql, = 7. We then fit for 
the linear coefficient, a^, for each of the cosine functions, 
so that the fitted model is: 

^ We thank Ron Gilliland for useful advise on this topic 



N 

X(t,) = ^aJ,(i,) (2) 

1=0 

We then subtracted model M from the lightcurve (in- 
cluding the transits). The model is shown over the data 
for the QO and Ql photometry in Figure [TJ 

2.3. Median Normalizations 

A second stage of normalization is applied to the 
data after the long-term detrending. Here, we split the 
lightcurve up into 18 individual transit and occultation 
events (giving 36 arrays in total) . Each array spans from 
-0.125 to -1-0.125 in orbital phase surrounding the event 
in question. The fluxes and associated errors in each ar- 
ray and then divided by the median o f each array. This 
is similar to the technique adopted bv lKipping fc BakosI 
(.20101 

2.4. Outliers 

Despite the PA processing, some outliers still remain in 
our detrended, normalized photometry. We must remove 
these before it is possible to perform the final lightcurve 
fits. Since these outliers can occur within the transit 
event itself, it is necessary to perform a preliminary fit of 
the transits and then remove outliers from the residuals. 

For the purpose of the identifying outliers, we perform 
an initial global fit (as described later in § 13.21) . The 
residuals are then used to search for outlier points by 
flagging those which occur 3-(t away from the model. 

2.5. Time Stamps 

In the DR5 handbook, the following advise is given: 
"The advice of the DAWG [Data Analysis Working 
Group] is not to consider as scientifically significant rel- 
ative timing variations less than the read time (0.5s) or 
absolute timing accuracy better than one frame time 
(6.5s) until such time as the stability and accuracy of 
time stamps can be documented to near the theoretical 
limit." 

Relative time differences correspond to, for example, 
performing TTV (transit timing variations) and TDV 
(transit duration variations) on the Kepler data alone. 
Absolute time differences corresponds to, for example, 
performing TTV and TDV on the Kepler data plus all 
previously observed data. We stress these limitations 
early on in our study. 

The Kepler time stamps from DR5 are in BJDutc 
(Barycentric Julian Date in Coordinated Universal 
Time) and a correction to BJDtdb (Barycentric Ju- 
lian Date in Te r restri al Dynamic Time) is advocated by 
lEastman et al.l ()2010l) in all transit timing studies. The 
correction between UTC and TDB is given by BJDtdb = 
BJDutc + N + 32.184, where A^ is the number of leap 
seconds which have elapsed since 1961. This correction 
has been applied all data analyzed in this study. 

2.6. Correlated Noise 

Time-correlated noi se may affect the estimation of 
lightcurve parameters (jCarter et al.ll2010l ) and so we here 
discuss the degree to which this data set is affected by 
correlated noise. We present two methods of assessing 
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th e degree of red nois e, following the approach adopted 
bv lCarter et"aI1 (|2009l ). 

First, using the residuals of our final fits (which will be 
introduced in more detail in i|3.2p . we bin the residuals 
into a bin size N and evaluate the r.m.s. of the data. We 
repeat this process from iV = 1 up to = 360 (which 
is approximately equal to the time span of each discrete 
lightcurve array, 0.25 days) and the results are shown in 
Figure O The figure reveals that our corrected data is 
follows closely the expectation of independent random 
numbers, ctat = aiN^'^^^[M/{M - 1)]^/^, where M is 
the number of bins. 

Second, we computed the lAUaiil (|19660 variance (t\{1) 
of the residuals, defined as: 

= 2(A^ + l-20 g [l^J'^' ' "'^'^j 

(3) 

where denotes the residual of the k^^ data points, 
N is the number of data points and / is the lag. For 
independent residuals, one expects ~ a\{0)/l, for 

which our residuals can be seen to be satisfy in Figure [21 

The bulk r.m.s. of our data set is 237.2 ppm. In com- 
parison, the PCD corrected photometry has a bulk r.m.s. 
of 230.5 ppm (after removing outliers). It is possible that 
PCD overfitted the data or that our own correction is an 
underfit. Based upon the analysis above though, we find 
no strong evidence for correlated noise in our corrected 
photometry, which would be expected for underfitted de- 
trending. 

3. MODEL DETAILS 

3.1. Model Generation 

3.1.1. Lightcurves 

The primary transit l ightcu rve model is computed us- 
ing the lMandel fc Agol (|2002[ ) limb darkening algorithm. 
The outputted fluxes are corrected for variable baseline 
flux, OOT, to propagate the baseline r.m.s. into the 
lightcurve. The occultation lightcurve is computed in 
the same way, except the limb darkening coefficients are 
fixed to zero and the final lightcurve is then squashed 
by a factor which is equal to the ratio of the transit to 



occultation depth. We note that multiplying the ratio- 
of-radii squar ed, jp^, by this factor a nd then feeding this 
value into the lMandel fc Ago! (2002*) code instead would 
be erroneous, since the algorithm would think the planet 
was very small leading to sharper ingress/egress features. 
By applying the transformation at the end, we preserve 
the correct lightcurve morphology. 

The true anomaly is calculated from the time stamps 
by solving Kepler's Equation at every instance. Tran- 
sit dura t ions are computed using the expressions of 
iKippind ()2008[ ). which account for orbital eccentricity. 
Although TrES-2b is be lieved to be on a circular orbit 
(jO'Donovan et al.ll2007D . using the most general equa- 
tions allows us to float the eccentricity parameters to 
propagate their uncer tainties. Recent Spitzer occulta- 
tion measurements bv^ O'Donovan et al.l (|2010D strongly 
constrain ecosw = (0.00053 ± 0.00102). We allow both 
esinoj and ecosw to be fitted for in our global fits, but 
as the e cos uj term move s away from the value found by 
lO'Donovan et al.l (|20100 . a penalty is assigned (see 
Equation [71) 

We initially use the errors from the normalized PA 
lightcurve and then rescale the errors after one iteration 
of the global fit. Errors are scaled such that the func- 
tion for the transit data, the occultation data, the RV 
data and the transit times are each equal to the number 
of data points in that fit minus the number of degrees of 
freedom in the model. 

3.1.2. Radial velocity 

The radial velocity curve is computed assuming a single 
planet in a Keplerian orbit. The free parameters in the 
model are the time of transit, the orbital period, ecosw, 
esinw and the semi-amplitude K. We do not consider 
the Rossiter-McLaughlin (RM) effect since our principal 
goal is to characterize the orbit and the points for the 
RM lead to very little improv ement in the parameters 
listed here (|Winn et al.ll2008D . but severe increases m 
CPU time. As a result , we o nly use the radial velocities 
from lO'Donovan et al.l (|2007( ). taking care to convert the 
times to BJDtdb- 

3.1.3. Transit times 

So far we have three data sets which are fitted for; 
the transits, the occultations and the radial velocities. 
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Bin size Lag 

Fig. 2. — Assessment of correlated noise. Left panel: The r.m.s. of the time-binned residuals as function of bin-size. Right panel: The 
square-root of the Allan variance of the residuals, as a function of lag. The gray lines are computed from the data and the black lines show 
the expected behavior for uncorrelated noise. 



Usefully, TrES-2b is a relatively old discovery and several 
years of transit measurements exist. However, most of 
these come from amateur measurements, which have not 
been peer-reviewed, and thus may not be as reliable. A 
resolution to this is to use median statistics to define the 
merit function and therefore provide a robust estimation 
of the goodness of fit, even in the presence of outliers. 
We therefore add in a fourth data set to our global fits 
coming from the timings, which provides extremely tight 
constraints on the ephemeris. 

Let us consider the typical merit of function first, which 
is based on mean statistics. In order to compute the 
ephemeris, we throw in a trial model of t + nP and then 
calculate the residuals for each point, r^. We then evalu- 
ate the weighted squares of each of these measurements, 
given by (ri/oiY , where Oi is the measurement error. In 
a normal analysis, we would then sum these weighted 
squares together to give the and then perturb the 
model until we obtain the lowest possible y^: 



i=l 

X = n — — 



= ?iMean{(r,/cr,)^} 



(4) 



Inspection of the above equation reveals the simple way 
in which we can change our merit function to be in the 
form of median statistics, to give : 



erature and the ETD (Exoplanet Transit Database^) are 
almost always in HJDutc (Heliocentric Julian Date in 
Coordinated Universal Time). We use the JPL Hori- 
zons ephemeris to convert the HJDutc times to BJDutc 
and then apply the correction for leap-seconds to yield 
BJDtdb- The list of used transit times is presented in 
the appendix, Table O 



3.2. Fitting Algorithm 

Fits are accomplished by using a Metropolis- 
Hastings Markov Cha i n Monte Carlo (MCM C) algorithm 
(jTegmark et al.l 120041 : iHolman et al.l 120061 ). The routine 
begins from a starting point, which we select to be 5- 
(7 away from the estimated solution, and then generates 
new trial parameters by making a jump computed us- 
ing a Gaussian proposal distribution centered upon the 
current position with a standard deviation given by the 
"jump size". Jump sizes are selected, usually through a 
process of iteration, to be equal to the 1-a uncertainties 
for each parameter. 

The trial parameters are then used to produce a model, 
which is compared to the observations to produce the 
goodness-of-fit merit function, x^- Trials producing a 
lower than the current position are always accepted 
and the trial position becomes the current position, con- 
stituting an accepted jump. Trials producing a higher 
are accepted with a probability: 



= 71 Median{(ri /cTj)^} 



(5) 



The distribution is very similar to that of the x^ dis- 
tribution, but a scaling factor is required to make them 
equivalent. This factor is frequently required when con- 
verting median statistics to mean statistics; for example 
the standard deviation is given by 1.4286 times by the 
median-absolute-deviation and the error on the sampling 
median is 1.253 times the error on the sampling mean. 
In this case, the factor was computed using Monte Carlo 
simulations where we found the factor 2.26 provides the 
correct scaling. 

Times of transit minimum found in the exoplanet lit- 



P(accept) = exp(-Ax^/2) 



(6) 



where Ax^ is the difference in between the current 
position and the trial position. The algorithm stops when 
125,000 trials have been accepted and the first 25,000 
(20%) are discarded as burn-in leaving 10^ points for the 
posterior distributions . Our algori thm follows the same 
procedure detailed in iFordI (|2005f) . The overall merit 
function (see for details) is given by: 



http://var2.astro.cz/ETD/ 
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(7) 



We fit using 14 free parameters {t, P, p'^,Ti,b, e cos w, 
esinw, OOT, OOS, wi, W2, Fp/F^, K, 7}, which we 
elaborate on here, r is the time of transit minimum for 
the optimum epoch (that epoch which produces the min- 
imum correlation to P) and is defined as the instance 
when the planet-star sky-projected separation is mini- 
mized (note, this is frequently given the misnomer "mid- 
transit time"). P is the orbital period, is the ratio- 
of-radii squared and b is the impact parameter (defined 
as the planet-star sky-projected separation in units of 
the stellar radius at the instance of inferior conjunction). 

T is the transit duration between the instance of the 
planet's center crossing the stellar limb to exiting under 
the same condition. Ti is the "one-term" approximate 
expression for this parameter, given by Equation 15 in 
iKippingj (^2010a^ ( an exact a nalytic form for T is not 
possible, see iKippinel (|2010aD for details). Stellar limb 
darkening is accounted for using a quadratic limb dark- 
ening model, modeling the specific intensity as a function 
of /i: 



^1 



(8) 



where n is the cosine of the angle between the observer 
and the normal to the stellar surface, ui and U2 are 
known to be highly correlated in typical lightcurve fits 
(|Palll2008l ) and instead we opt to use wi and 'W2, which 
are related to the quadratic limb darkening coefficients 
(|Pa]|l2008n via: 



Wi = Ui cos tp — U2 sm tf 

W2 — U2 cos ip -\- Ui&inip (9) 

m (|2008f ) has advocated using this linear combination 
instead of ui and U2 due to the improved decorrelations. 
The author also recommends using Lp — 40°, which was 
done so in this study. During the MCMC, we discard 
any trials which yield unphysical limb darkening coeffi- 
cients, defined as those which are not everywhere positive 
and produce a monotonically decreasing profile from limb 
to center. This is implemented by using the conditions 
(jCarter et al.l l2009l: U1 + U2 < 1, U1+U2 > and ui > 0. 

Finally, the final parameter values quoted in this pa- 
per are given by the median of all of the accepted MCMC 
trials for the parameter in question. Similarly, 1-a uncer- 
tainties are calculated by evaluating the 34.1% quantiles 
either side of the median. 

3.2.1. Why fit for eccentricity? 

Some readers may question why we choose to fit for 
eccentrici ty when the orbit is c onsistent with a circu- 
lar orbit (jO'Donovan et aLll2007l) . Firstly, we point out 



that by using all of the known transit times, the Kepler 
lightcurves and occultation constraints we are able to de- 
rive the most precise constraints on e yet for this system, 
which is a worthwhile goal in itself. 

However, the most important reason for fitting for e is 
that any uncertainty on e leads to inflated uncertainties 
on the derived stellar density, p*. As pointed out by 
iKippind (|2010aD . the retrieved stellar density is given by 
the approximation ~ p,^circ/^(e, t^j) where the first 
term is the stellar density derived from a circular fit and 
5* is given by: 



* = 



(1 -I- e sinw)" 
(1 + e2)3/2 



(10) 



P*.circ is determined purely photometrically and thus 
the uncertainty will decrease as ~ l/VA'transits, where 
-^transits is the uumbcr of observed transits by Kepler. 
This parameter can therefore be expected to be known 
to very high precision by the end of the Kepler Mission, 
given the short orbital period of TrES-2b. In contrast, 
5" can only be measured by radial velocities and/or oc- 
cultation events. Given the visible bandpass of Kepler, 
occultation events are not expected to be detectable for 
the majority of transiting planets, and so the radial ve- 
locity determination dominates. With typical transiting 
planets receiving sparse radial velocity coverage, it can 
be appreciated that the uncertainty on 'J will often be 
the limiting factor in the measurement of a precise p* . 

The point is that we do not know the orbit is exactly 
circular (indeed this is practically impossible) and thus 
we cannot assume esinw = and ecosw = exactly. 
In reality, we have errors on both of these and can only 
say it is circular to within a certain confidence level. This 
uncertainty therefore propagates into a much larger error 
for th e stellar density. As an example, iKipping fc Bakoi 
compare fits for Kepler-4b through 8b using both 
circular and eccentric fits and find the errors on p» con- 
sistently infiate for the latter. 

3.2.2. Why fit for limb darkening? 

Another methodology we adopt, which is not a com- 
pletely standard practice in the exoplanet literature, is 
that we fit for the limb darkening coefficients. Fitting 
for quadratic limb darkening requires a very high signal- 
to-noise if one wishes to achieve convergence, especially 
for a near-grazing transit. In many ground-based mea- 
surements, it is not possible to fit for these coefficients, 
although linear limb darkening could be used instead. 

However, if fitting for the limb darkening is viable, it 
is always preferable. This is because transit parameters 
derived using fixed limb darkening coefficients are funda- 
mentally model dependent, where the model is that of the 
stellar atmosphere model. In contrast, transit parame- 
ters derived using fitted limb darkening are independent 
of a stellar atmosphere. This makes them vastly more 
robust and reliable. 

This point is particularly salient for TrES-2b. For a 
near-grazing transit, the planet only ever crosses the 
limb, where the star is most severely darkened. Thus 
the choice of limb darkening coefficients has a very sig- 
nificant effect on the derived planetary radius and transit 
depth especially. The total stellar fiux, which defines the 
observed transit depth, is essentially extrapolated from 
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the stellar centre to the limb based upon the fitted limb 
darkening coefficients of the limb only. This leads to large 
correlations between the limb darkening coefficients and 
the ratio-of-radii squared. 

3.3. Blending 

Recently, iDaemgen et al.l ()2009[ ) showed that the 
TrES-2 has a very nearby star, which was proposed 
to be in binary star system composed of the originally 
known GO TrES-2A star and a previously undetected 
K4.5-K6 companion, (labeled TrES2/C by the authors). 
In the z'-band, the magnitude difference was estimated 
to be 3.43 and thus we estimate the blending factor 
B (which is defined in iKipping fc Tinettil (|2010h ) to be 
B= (1.04246 ± 0.00023). 

This blending acts to dilute the transit depth and thus 
causes us to underestimate the true planetary radius. 
Correcting for blends may be accomplished by follow- 
ing the prescription of IKipping fc Tinettil ()2010f ). which 
we adhere to in this work. Self-blending due to nightside 
emission is expected to be negligible in the Kepler band- 
pass (see same work) and thus need not be accounted 
for. 

3.4. Limb Darkening Computation 

In §4.11 we will discuss how limb darkening coefficients 
are fitted for in the final results. However, it is useful 
to generate the limb darkening coefficients from theoret- 
ical models for i) providing a sensible starting point for 
the fitting procedure ii) later comparison of theoretical 
models versus fitted limb darkening. 

Limb darkening coefficients were calculated 
for the Kepler bandpass for TrES-2b. For 
the Kepler bandpass, we used the high reso- 
lution Kepler transmission function found at 
[htt p : / /keplergo . arc . nasa. gov / CalibrationResponse . shtml , 
We adopted the SME-de rived stellar properties r eported 
in ISozzetti et al.l (|2007f ). We employed the iKuruczl 
(|2006l ) atmosphere model database providing intensities 
at 17 emergent angles, which we interpolated linearly 
at the adopted T^s and log 5 values. The passband- 
convolved intensities at each of the e mergent angles were 
calculated following the procedure in IClaretl (|2000[ ). To 
compute the coefficients we used the limb darkening law 
given in Equation [S] 

The final coefficients resulted from a least squares sin- 
gular value decomposition fit to 11 of the 17 available 
emergent angles. The reason to eliminate 6 of the angles 
is avoiding excessive weight on the stellar limb by us- 
ing a uniform sampling (1 fj, va lues fro m 0.1 to 1, plus 
H = 0.05), as suggested bv lDfaz- Cordov cs et al.l (fl99^) . 

3.5. Drifts and Trojans 

Before we provide the final results, we discuss how we 
performed global fits including a linear drift in the radial 
velocities, 7, and a temporal offset between the radial 
velocity null and the time of transit minimum. At (such 
a temporal offset is e xpected to be induced by Trojans 
(|Ford fc Gaudill2006l )). By switching on and off these 
parameters, there are four possible permutations of the 
fits we can execute; eight when one switches on/off ec- 
centricity as well (see Table [T]). 

In general, fitting for an excessive number of free pa- 
rameters is undesirable as it increases the errors on the 



TABLE 1 

Bayesian Information Criterion (BIC) values for eight models 
executed for the global fits. We highlight the lowest BIC values. 



Model BIC 



e = 0, 


7I 


= 0, 


Ai| 


= 


34048.7 


e = 0, 


7I 


> 0, 


Af 


= 


34064.2 


e = 0, 


7I 


= 0, 


Ai 


> 


34058.6 


e = 0, 


7I 


> 0, 


At| 


> 


34069.7 


e > 0, 


7I 


= 0, 


M\ 


= 


34067.2 


e > 0, 


7I 


>o, 


At 


= 


34077.5 


e > 0, 


7I 


= 0, 


Ai 


> 


34259.1 


e > 0, 


7I 


> 0, 


At 


> 


34268.4 



other terms. In order to decide whether these two 
additional parameters should be included or not, one 
may evaluate the Bayesian Inform ation Criterion (BIC) 
(|Schwarzlll978l: iLiddie et al]|2007D . for each of the pro- 
posed models. The model with the lowest BIC is ac- 
cepted and subsequently used in the global fits reported 
in the next section. These fits used 125,000 MCMC trials 
with a more aggressive minimization downhill simplex 
implemented afterwards. It is based on this lowest so- 
lution from which the BIC is computed. 

We therefore performed eight versions of our global fits, 
with the results for the BIC values presented in Table [TJ 
We find that neither a drift nor a temporal offset are 
accepted for either the circular or eccentric models. We 
therefore proceed to consider them fixed to zero. The 
results, however, do allow us to place upper limits on 
7 and At. We find I7I < 0.12 ms~ Vear~^ and |At| < 
0.15 days to 3-(t confidence. This excludes a > 0.94 Mj 
Trojan in a 1:1 resonance with TrES-2b. 

4. RESULTS OF GLOBAL FITS 

The global fits were performed using the full Kepler 
time series as described in §3.11 The final results are 
given in Table [H After the main MCMC fits, a down- 
hill simplex routine is used to obtain the lowest solu- 
tion. We plot this solution over the data in Figure [3f. 
Histograms of the marginalized posterior distributions 
for each of the fitted parameters are shown in Figure |4l 
which clearly indicate convergence of the fitting param- 
eters. 

4.1. Limb Darkening Fitting 

Fitting for the limb darkening (LD) coefficients is chal- 
lenging because TrES-2b is a near-grazing transit and 
thus only samples a fraction of the stellar surface. How- 
ever, the extremely high quality of the Kepler SC pho- 
tometry and the fact we have 18 transits does allow 
for a good solution (given in Table [2]) . The inevitably 
strong correlations between the quadratic coefficients is 
presented in Figure [5] 

We find that the theoretical limb darkening coefficients 
lie within the 1-cr confidence regio n of our fi ts, i ndicating 
an impressive prediction for the iKurucd (|2006i) atmo- 
sphere model. One major benefit of fitting for the limb 
darkening is that the uncertainty in the stellar properties 
is built into the model and thus leads to larger, and ulti- 
mately more realistic, estimates of the various parameter 
uncertainties. Parameters which are highly correlated to 

^ A high definition version of this figure is available at 
www.homepages.ucl.ac.uk/~ucapdki/globalfit.pdf 
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TABLE 2 

Results from global fits of TrES-2b using eighteen short- cadence (SC) Kepler transits. We show results for both circular and eccentric fits 

in columns 2 and 3. In column 4, we provide the previous estimates of the system parameters from \0' Donovan et al\ I200'\f 
\Holman et a l. {2001)" and \Sozzetti et oil i200^Y " . In general, the eccentric fit leads to more realistic errors. Quoted values are medians 
of MCMC trials with errors given by 1-a quantiles. * = fixed parameter; f = parameter was floated but not fitted. 



Parameter 



Circular 



Eccentric 



Previous 



Model indep. params. 



^■4 'Ubi«yb_o. 00000016 



0.00000012 
t-0. 000022 
-0.000021 
25 



P [days] 
[BJDtdb - 2,450,000] 
Ti,4 [s] 
fi [s] 
72,3 [s] 

(Ti.a ~ T3,4) [s] 
(Rp/R* f [%] 
b 

^occultation [PPm] 

esinoj 
ecosoj 

K [ms-i] 
7 [ms"-*-] 

B 

ui 

U2 

Rp/R* 
a/R* 



P* [g/,cm ^] 
log(9p [cgs]) 



4849.526635 
6438 
4624 



0.000026 
0.000026 
+31 
33 
+42 

1950tlT« 
2242lg 
1.643l°;°i 

91+23 
^^-22 

0* 
0* 
1* 

181.4l^-f 
-29.9lfj 
1.04246 ± 0.00023 f 



2.47061892 
4849.526640 



6439_2g 
4624ti^j 
1942t«4 

2247^37 

1 CQQ+0.076 
1.000_Q 

O.S41S_Q 0045 
19+18 



-0.009 
0.0005 
0.973 
181.0; 
-29.2 



17 

+0.024 
0.029 
+0.0018 
0.0018 
+0.071 
-0.082 
,+5.5 
5.4 
+2.6 
2.6 



0.52 



+0.44 



0.34 
+0.37 
-0.48 
+0.0032 
-0.0020 
7 QH'^+'^ -'^'^^ 
I .»oo_o 084 

83.952tO;J3i 



0.06 
0.1282 



1.04246 ± 0.00023 f 
45+0-42 
U.^O_0.30 

r, + 0.33 

-0.46 



0.12"' 



0.1278 
8.06 



+0.0029 
-0.0018 



0* 



1.63"' 



,+0.16 
^-0.13 
q Q1 7+0.018 
( _o,oi8 



+0.25 
-0.21 

84.07lg:^t 

o.oi8t°:«f3 

2681180 

^■°'^-0.13 
q n9f.+0.025 
o.ozo_o 024 



2.470621 ± 0.000017 
6624 ± 72 " 

2459 ± 162 " 

0.8540 ± 0.0062 " 

0* " 
0* 

1* " 
181.3 ± 2.6 ' 

1* " 
0.22* " 
0.32* " 
0.1253 ±0.0010 " 
7.63 ±0.12 " 
83.57 ±0.14 " 
0* " 

1.375 ±0.065 



Model depend, params. 



Teff [K] (SME) 
Iog(g[cgs]) (SME) 
(Fe/H) [dex] (SME) 
M. [Mq] 

R* [R@] 

Iog(9 [cgs]) 

L, [Lq] 
My [mag] 
Age [Gyr] 
Distance [pc] 
Mp [Mj] 
Rp [Rj] 
pp [gcm-3] 
a [AU] 



5850 ± 50 t 
4.4 ±0.1 t 
-0.15 ±0.10 t 

u.yy^_o.o5o 
n Qc;s+0-0i8 

0.y5»_o.020 

4 469+0-015 
'i.'iO3_0 012 

U.yDi_o 057 
, „„-+0.074 
' -0.072 



3.3" 



-1.9 



-1.3 
202.6lg8 
1 r,oc+0.d58 
l.zuo_o 058 
1 199+0 020 
l.l»3_0 022 

O.STOiHS 



5850 ± 50 t 
4.4 ±0.1 t 
-0.15 ±0.10 t 
qqn+0 041 
U.M3U_o 048 

u.aoz_o 029 
4 475+0-024 
^•^'■-'-0.024 

948+°-°^^ 
u.a'io_o 070 

4 892+°-'""' 
'i.o»z_o 088 

Q 1+2.0 

•^•^-1.6 
201. 2^^-^ 

1 9^9+0-050 
l.iUZ_o 051 
1 1S7+0-034 
J^-^8<-_o.o35 
n SQ1 +0-074 
0.891_o 064 

U.UOODO_o 00058 



5850 ± 50 
4.4 ±0.1 
-0.15 ±0.10 
0.980 ± 0.062 

i.oooj 

4.426^ 



^+0.036 
-0.033 

r. + 0.021 

" -0.023 



4.77 ±0.09 

5 1+2.7 iii 

— 2.3 
220 ± 10 
1.198 ±0.053 
1.222 ± 0.038 

0.0367t°-°°i,? 



the limb darkening coefficients, such as the the transit 
depth (see §3.2.2p . have their associated errors increase 
considerably as a result of this process. 

As discussed in we actually fitted for wi and W2 
rather than ui and U2 to decrease the correlations, fol- 
lowing the prescription of lPall ()2008l ). We chose (p = 40°, 
as this was su ggested as a useful first-guess for the term 
by iPall ()2008l ). However, future studies of this system 
would benefit by using a more optimized value of ip. By 
using a principal component analysis (PCA), we are able 
to find this optimum angle to be i^j = 42.7033°, very close 
to the 35° - 40° range advocated by Pal (2008). 

It is important to consider the effects of fitting for LD 
carefully. We re-ran our fits with the LD parameters 
fixed to their best- value and found that the errors on nu- 
merous parameters were considerably reduced, in many 



cases by an order-of-magnitude. As an example, the 
transit depth error is reduced by a factor of 17.5 when we 
fixed the LD parameters. The errors found using fitted 
LD correspond to the absolute uncertainty in each pa- 
rameter. Therefore, if we wish to compare the duration 
found from Kepler photometry with, say, a ground-based 
measurement in a different bandpass, we must fit for the 
LD parameters separately in both cases. However, if we 
consistently employ the same bandpass and instrument 
response function for the same star, then there is no need 
to refit the LD parameters everytime. By fixing the LD 
parameters to their best-value, we compute the relative 
duration changes, within that bandpass. 

For TTV, the error in the time of transit minimum does 
not appreciably change between fitting and not-fitting 
the LD parameters. Therefore, the TTV seems to be re- 
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Fig. 3. — Short-cadence folded transit lightcurve of TrES-2b 
( circles ) with model fit overlaid. Residuals from the fit are shown 
below, offset by +0.982. 

liable across different bandpasses and instruments. This 
opportunity will be exploited later in [JT) 

4.2. OccuUation 

In the short cadence global fits, we do not detect an 
occultation for the planet in the circular nor the eccen- 
tric fits. Although no occultation is detected, a robust 
upper limit is obtained. We choose to use the eccentric 
fit from here on, as it provides the most realistic errors 
(see g3XT|). 

The best-fitted occultation depth is ^occultation = 
19^17 ppm, indicating no detected signal. The posterior 
distribution of the occultation depth is presented in Fig- 
ure SI We exclude an occultation of depth > 72.9 ppm 
to 3-CT confi dence. 

Recently, ISpie gel fc Burrow ^ (201(f) predicted that the 
occultation of TrES-2b, in Kepler's bandpass, would be 
< 20 ppm, assuming no refiected light contribution. Our 
results are therefore highly consistent with the theoreti- 
cal models for this planet. 

The 3-<7 limit constrains the geometric albedo to be 
Ag < 0.146 and a dayside brightness temperature of 
Tp^day < 2413 K (for comparison, the equilibrium tem- 
perature is 1472 K). We note that our 3-ct lin iit is tighter 
than t hat for HD 209458b as measured by iRowe et al.l 
(|2008') using MOST, where Ag < 0.17 to 3-cr confidence. 
Therefore, TrES-2b is currently the darkest exoplanet 
known to exist. For comparison, the upper limit corre- 
sponds to a planet of similar albedo to Mercury (0.138). 

4.3. A Search for Asymmetry 

Lightcurve asymmetry is generally not expected but 
may reveal interesting, new physics for hot- Jupiters. One 
possible source would be an oblate star with a signifi- 
cant spin-orbit misalignment causing an asymmetry in 
the ingress/egress durations. We here describe how we 
searched for asymmetry in the lightcurve. 

We divide the folded lightcurve into points before and 
after the globally fitted time of transit minimum, where 
the fold is performed using the globally fitted period. We 
then mirror the two halves upon each other to search for 
signs of asymmetry in the lightcurve. The residuals of 
these two halves are shown in Figure [6] 



We first perform a linear interpolation of the folded 
data prior to the time of transit minimum. This func- 
tion is then evaluated at the time stamps of the folded 
data after the time of transit minimum and the associ- 
ated uncertainties are carried over. We then subtract the 
two and add the two sets of fiux uncertainties in quadra- 
ture. The "mirror residue" exhibits an r.m.s. scatter of 
305 ppm, whereas from a theoretical point-of-view one 
expects scatter equal to x 237.2 ppm = 335 ppm. A 
chi squared test gives 6990 for 7684 data points. The 
ingress and egress therefore exhibit remarkable symme- 
try. 

The most significant feature in Figure IH] is a slight drop 
at around +0.02 days. This feature is only 2-a significant 
with the current data, but could be scrutinized further 
in later data releases. 

4.4. Eccentricity 

As shown in Table [51 we performed both circular and 
eccentric fits to illustrate the consequences of fitting for 
eccentricity. The fits find very similar values, with 
the circular fit being marginally larger by Ax^ — 1.89 
for 30697 data points. Using an F-test, we find that the 
eccentric fit is accepted over the circular fit with a con- 
fidence of 38.9%, which we consider to be insignificant. 
Therefore, we conclude the orbit of TrES-2b is consistent 
with a circular orbit, based upon the current data. Fur- 
ther, using the marginalized posterior distribution, we 
estimate that the eccentricity satisfies e < 0.094 to 3-a 
confidence. 

The e cos w prior from lO'Donovan et all (|2010D places 
a much stronger constraint than that obtained for either 
component purely from the radial velocity. As a result, 
we find a much larger uncertainty on e sin w than e cosoj. 
Whilst both components are consistent with zero, it is 
unlikely from an a-priori perspective than e sin uj will 
be non-zero given that ecosLJ is essentially zero. If 
both components had similar uncertainties, but consis- 
tent with zero, then the eccentricity would be tied down 
to e < 0.0085, but we stress that this is not a conclusion 
which can supported purely based upon the data. 

4.5. Revised Masses and Radii 

Fundamental parameters of the host star such as the 
mass (M,) and radius (i?*), which are needed to infer 
the planetary properties, depend strongly on other stel- 
lar quantities that can be derived spectrosco pically. For 
this w e used the spectroscopic analysis of Sozzetti et al.l 
(|200l who determine Toff = (5850 ± 50) K, [Fe/H] = 
(-0.15 ± 0.10) and log(g [cgs]) = (4.4 ± 0.1). 

In principle the effective temperature and metallicity, 
along with the surface gravity taken as a luminosity in- 
dicator, could be used as constraints to infer the stellar 
mass and radius by comparison with stellar evolution 
models. 

For planetary transits a stronger constraint is often 
provided by the a/i?* normalized semi-major axis, which 
is closely related to p*, the mean stellar density. The 
quantity a/R^, can be derived di rectly from the transit 
lightcurve (Se ager fc Mallen-Orne las 20 03) a nd the RV 
data (for eccentric cases, see iKippingl (|2010al )). The re- 
sults of our 100,000 MCMC trials are used to produce 
an array of 100,000 estimates for p*. Toff and [Fe/H]. For 
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Fig. 4. — Marginalized parameter posterior distributions. Reading down from top left to right, we have t, P, , T\, b, OOT, i5occultation) 
OOS, wi, W2, K, 7, esint*; and ecosw. Results come from global fits assuming no Trojan body and no linear drift in the RVs, but allowing 
orbital eccentricity. 



every tri a l, we match stellar evolution isochrones from 
lYi et alj (|2001[ ) to the observed properties to produce 
100,000 estimates of the absolute dimensions of the star. 
Finally, the planetary parameters and their uncertainties 
were derived by the direct combination of the posterior 
distributions of the lightcurve, RV and stellar parame- 
ters. 

After the first iteration for determining the stellar 
properties, as described in iBakos et al., (,2009, ) . we find 



that the surface gravity, log((ji [cgs]) = 4.475^q q24' is 
highly consistent with thelS ozzetti et all ()2007l ) analysis. 
Therefore, a second iteration (which would use the new 
log (7 value) of the isochrones was not required and we 
adopted the values stated above as the final atmospheric 
properties of the star (shown in Table [2]). 

The revised parameters are in excellent agree - 
ment with the estima tes fro m lO'Donovan et al.l (l20Q7t ). 
iHolman et al] (|2007[ ) and ISozzetti et all (|2007D . all 
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Fig. 5. — Distribution of the quadratic limb darkening coefficients 
from the MCMC global fits of the short-cadence data. Black points 
correspond to the 3-cr region and gray to the 1-a. The diamond 
marks the the oretical limb darkening coefficients computed from a 
\Kuruci (20061) style atmosphere. 
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Fig. 6. — Mirror residuals. Subtracting the two halves of 
the lightcurve from each other produces the "mirror" residuals, 
which would reveal any signs of lightcurve asymmetry - how- 
ever, none are evident. The vertical black lines mark the start 
of ingress/egress and the end of ingress/egress. 

shown in Table [2] for comparison. However, our de- 
rived stellar density is markedly larger and this leads to 
a slightly smaller, more massive star, which consequently 
'deflates' TrES-2b slightly. 

5. TRANSIT TIMING VARIATIONS (TTV) 

We will here only consider short-term transit timing 
variations, which we define to be those occurring within 
the timescales of the eighteen observed Kepler transits 
in QO and Ql. A long-term transit timing analysis is 
provided in SJT) 

5.1. Fitting Method 

For the individual fits, we do not expect limb darkening 
to vary from transit to transit and thus using a single, 
common set of LD coefficients is justified (as explained 
earlier in §??). We therefore fix the quadratic coefficients 
to those found to give the lowest in the global fit 
we performed earlier (selected values were wi = 0.22366 
and 'W2 — 0.39288). Aside from the limb darkening, the 
eccentricity terms e sinw and e cos w are also held fixed to 



the lowest solution (esinw = —0.020864 and ecosw = 
0.00076088). In total, there are five free parameters used: 
{r, p^, Ti, 6, OOT}. An initial run is used to compute 
scaling factors individually for each transit epoch, which 
span ±0.125P of the linear ephemeris predictions. The 
scaling factors are selected so that the lowest solution 
found is equal to the number of degrees of freedom in the 
model, as before. 

The individually fitted transit lightcurves are shown in 
Figure [7] and the parameters in Table [31 We note that 
in none of the transits is a second transit-like feature 
observed, as claimed bv lRaetz et al.l (|2009D . 

5.2. Control Data 

We describe here how we produce control data in the 
form of artificial lightcurves. The act of producing con- 
trol data by which to compare the genuine observations 
is a practice frequently applied in many aspects of sci- 
entific study. In our case, the control data serves two 
principal functions: 

1. Rescaling of the parameter uncertainties 

2. Identification of signals due to "phasing" 

5.2.1. Rescaling 

Both of the s e iss ues were first noted in 
iKipping fc Bakod ()2010[ ). although control data was 
not used. The rescaling issue was observed by the 
authors as they found that the errors produced by the 
Metropolis-Hastings MCMC method led to scatters in 
their TTV and TDV much lower than the parameter un- 
certainties. The chance of this occurring by coincidence 
in an isolated case was estimated to be ~10%, however 
the pattern was recurring for the majority of parameters 
evaluated. This led the authors to conclude that there 
was strong evidence the measurement uncertainties were 
being overestimated. 

Calculating the necessary rescaling factor can be 
achieved by generating control data. For a single global 
fit, as performed earlier, this would be too time consum- 
ing with the 30,000-t- data points plus the correlations be- 
tween limb darkening and depths taking several weeks to 
fit in a single run. Therefore, the uncertainties presented 
in Table [5] may in fact be overestimates as well, although 
we have not confirmed this. For the individual data, one 
may take advantage of the fact that a planet exhibiting 
no TTV, TDV, T^V (depth variations), TbV (impact 
parameter variations) or baseline variations should yield 
a equal to the number of degrees of freedom in each 
case. For example, for the TDV, we have 18 transits 
with one model parameter, the mean duration, and so 
we expect = 17. 

Since no real system can be assumed to be absolutely 
temporally invariant, the only practical way forward is 
to generate artificial data for the control. To accomplish 
this, we take the global fit model found earlier and sam- 
ple it at the exact time stamps in each individual transit 
epoch array. The global model implicitly assumes that 
no parameters vary from epoch to epoch, satisfying our 
control condition. Next, we introduce Gaussian noise 
into the lightcurve equal to the actual noise recorded at 
those time stamps. This noise includes the scaling factors 
found in the individual fits, to ensure equivalence. These 
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Fig. 7. — Individual transits of TrES-2b from QO and Ql Kepler photometry. Top lightcurve is Kepler epoch going sequentially in 
time to epoch 18 at the bottom. Residuals are shown below. 
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control lightcurves are then fitting using the Metropolis- 
Hastings MCMC method using identical starting posi- 
tions, jump sizes, etc as the individual fits. 

The various parameter variations are then evaluated 
and the necessary scaling factor is computed. The scal- 
ing factors are given in the last line of Table |3l Our 
result s agree with the conclusions of iKipping fc Bakoi 
()201Cl() i.e. that all parameters have overestimated errors 
by around a factor of ~ 2. The reason for this overes- 
timation is unclear and despite close examination of our 
routines, we can find no obvious reason why this should 
occur. An indepen dent code used for HATNet discoveries 
(jBakos et al.ll2009t ) , which also uses Metropolis- Hastings 
MCMC, finds very similar uncertainties to the algorithm 
used in this work (detailed comparison s of the two meth- 
ods h ave b een previously provid ed in IKipping fc Bakoi 
(|201Cl[ ) and IKipping et al] ()2011|) ). suggesting this is not 
a specific flaw in our routine. 

5.2.2. Phasing 

In IKipping fc Bakod (|201Cl[ ). the authors considered 
a new term which they labeled as the transit "phas- 
ing". This corresponds to the time difference between 
the expected time of transit minimum and the nearest 
data point. For example, for data of cadence 60s, we 
would expect this time difference to be in the range 
±30 s. Phasing does not seem have a linear correla- 
tion to observed parameter variations, but does introduce 
false periods into the p ower spectrum of the variations 
(|Kipping fc Bako"il2010[ ). Removing the phasing effects 
is not currently possible, but we can at least generate the 
effects which phasing induce to compare to the real data. 

By generating our control data with the exact same 
cadence and time stamps as the data which we fit in 
the individual transit arrays, we can recover any possible 
influence the phasing may have on our results. In what 
follows, figures showing parameter variations always have 
the real data on the left-hand-side and the control data 
on the right, so the effects of phasing are most clearly 
visible. 

5.3. Analysis of Variance for TTV 

The TTV, shown in the top- left panel of Figure |H1 ex- 
hibits a r.m.s. scatter of 5.16 s, which demonstrates the 
impressive precision of these Kepler measurements. Af- 
ter rescaling the uncertainties, the scatter in the data is 
consistent with a linear ephemeris, exhibiting a — 19-0 
for 16 degrees of freedom. The excess scatter is l.l-cr sig- 
nificant, which we consider below our detection thresh- 
old. The unsealed errors yield — 6.2 for 16 degrees of 
freedom, supporting the hypothesis that the errors are 
significantly underestimated and justifying our rescaling 
methodology. Figure |8| shows the results. 

5.4. F-test Periodogram for TTV 

The F-test periodogram fits sinusoidal waveforms 
through the data of various periods, stepping through 
from the Nyquist frequency to the observational window 
in equally spaced steps of size 1/1000 of an epoch. Fit- 
ting for amplitude and phase, the is computed in each 
step, and then the F-test is performed. The false-alarm- 
probabilities (FAP) of these F-tests are then plotted in 
a periodogram. It is important to appreciate that the 



F-test is designed to look for sinusoidal waveforms, and 
thus periodic but non-sinusoidal waveforms would have 
their significances attenuated. 

The control data reveals periodogram peaks at 2, 4 and 
8 cycles which are harmonics of the sampling cadence of 
one transit measurement per transit epoch. In the real 
data, only one peak surpasses 95% confidence occurring 
at a period longer than the observation window. Such 
long period peaks cannot be considered genuine unless 
further transit epochs confirm the periodicity. In conclu- 
sion, there is no evidence for a TTV signal in the Kepler 
QO and Ql TrES-2b photometry 

5.5. Excluded TTV Signals 

We conclude our analysis of the TTV by evaluating the 
constraints on other planets, moons and Trojans in the 
system. For 16 degrees of freedom, r.m.s. scatter pro- 
ducing a = 36.2 is excluded to 3-cr confidence. This 
excludes r.m.s. scatter of 7.11s to the same confidence 
level. 

An outer perturbing planet in a j:j + 1 mean mo- 
tion resonance (MMR) would cause the inner transiting 
planet to librate lead ing to TTVs (jHolman fc Murray! 
120051 : 1 Agol et al.ll2005D . For 1:2, 1:3 and 1:4 resonances, 
the libration periods are 18.1, 10.5 and 7.2 cycles respec- 
tively. We therefore possess sufficient baseline to look for 
all such resonant planets. This excludes the presence of 
coplanar, MMR planets in these resonances of 0.11 M^, 
O.17Af0 and 0.22 Af® respectively. 

For an extrasolar moon in a retrograde orbit, the max- 
imum dYiiamkall5;_&^ separation is 0.9309 Hill 
radii (jPomingos et al.l [20061 ). For such a body on a cir- 
cular, coaligned orbit, we are able to exclude moons of 
1.15 A'/q. As the orbital separation decreases, we are able 
to exclude moons of masses > (1.07//) sini^ Mq, where 
/ is equal to the moon's orbital separation in units of 
Hill radii. 

Trojan bodies can also induce TTVs and thus con- 
straints on their presence can be also e stablished. Us- 
ing Equation 1 of Ford fc HolmanI (|2007( ) , and assuming 
a Trojan of angular displacement 10° from the La- 
grange point, we are sensitive to Trojans of cumulative 
mass > 0.46 to 3-cr confidence. However, the ex- 
pected libration period would be ~ 75 cycles and thus 
we do not yet possess sufficient baseline to definitively 
exclude such bodies. 

5.6. Proposed 0.21 Cycle Period Signal 

Another s ignal we are able to investigate is the one 
proposed bv lRabus et al.l (|2009l ). The authors claimed a 
0.21 cycle period sinusoid of amplitude 50 s provided a 
best-fit to the previously known transit times of TrES- 
2b, with a FAP of 1.1% and suggested a 52 exomoon 
as a possible origin. 

Fixing the amplitude and period to the proposed value 
and fitting for the phase term, we find a x^ = 1610 for 
the 18 data points. In contrast, the static model obtains 
a x^ = 19.0, which therefore excludes the claimed signal 
to high confidence. This highlights the dangers of looking 
for signals below the Nyquist frequency. 

6. TRANSIT DURATION VARIATIONS (TDV) 
6.1. Choosing a Statistic 
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TABLE 3 

Transit parameters of TrES-2b from individual fits of the SC data. Kepler epoch is defined as the first transit observed by Kepler. In 
this data, we fix the limb darkening to tfie best fit values from the global fit. Therefore, the data can be used to look for relative changes 
between these 18 tmns 'its, but not against previous ground-based measurements of TrES-2b. Errors have not been re-scaled, but the 
scaling factor is provided on the last line, based upon the analysis of control data. 
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Fig. 8.- - TTVs of TrES-2b. Top left: Observed TTVs for TrES-2b using the ephemeris of the global fit. Bottom left: F-test periodogram 
of the observed TTVs. Top right: TTVs computed from control data (artificial lightcurves) . Bottom right: F-test periodogram of the 
control TTVs. 
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Due to the near- grazing nature of the transit, the stan- 
dard assumption that T is th e optimum statistic for TDV 
searches may not be vahd (jCarter et al.ll2008t iKippingj 
[2003). In particular, the first-to-fourth contact duration, 
Ti^4, could potentially offer greater sensitivity. We found 
that the typical error on Ti_4 was 0.28% in the individual 

fits, whereas Ti marginally better at 0.24%. Another fac- 
tor in choosing a statistic comes from the effects of limb 
darkening. Whilst here we fix the limb darkening co- 
efficients to the best-fit values from the global MCMC, 
it is preferable to still avoid using a statistic which is 
strongly correlated to limb darkening. This is because 
we may be using slightly incorrect limb darkening coef- 
ficients which therefore feed into incorrect duration esti- 
mations. Whilst this is generally unavoidable, strongly 
correlated terms would clearly excerbate the situation. 
We find that Ti 4 has a correlation of 0.021 to the limb 
darkening coefficient W2 (the most strongly constrained 
coefficient) whereas Ti has a correlation coefficient of 
0.81. On this basis, we choose to use the Ti^4 statis- 
tic in what follows, defining the TDV of the i**^ transit 
measurement as: 

TDV, = [Ti,4], - [Ti,4]global (11) 

6.2. Analysis of Variance for TDV 

The TDV, shown in the top-left panel of Figure [HI 
exhibits a r.m.s. scatter of 22.4 s. After rescaling the 
uncertainties, the scatter in the data is consistent with 
a constant duration, exhibiting a = 15.4 for 17 de- 
grees of freedom. The unsealed errors yield — 5-6 for 
17 degrees of freedom, again supporting the hypothesis 
that the errors are significantly underestimated. Figure |9] 
shows the results. 

6.3. F-test Periodogram for TDV 

We continue by computing the F-test periodogram for 
the TDV data (shown in lower-left of Figure ^. The 
TDV data yields only one interesting peak occurring with 
a broad distribution surrounding 4.72±0.10cycles, signif- 
icance 94.1%. Firstly, this is below our formal detection 
threshold. Secondly, the peak seems to occur in the con- 
trol data, with distinct phasing periods occurring at 3, 
5 and 8 cycles. In light of this, we do not consider the 
signal to be genuine. 

6.4. Excluded TDV Signals 

The TDVs exclude a signal of r.m.s. amplitude 35.1s 
to 3-cr confidence, or variations in the duration of 0.77% 
over the 18 cycles. This excludes exomoons inducing 
TDV-V of mass > 23.5^/f cos is to the same con- 
fidence level. Additionally, it excludes fMssmis > 
1.17 Mi:^ through the TIP effect. 

Combining the TTV limits, the TDV-V limits and the 
TDV-TIP limits allows us to plot the parameter space 
of excluded exomoon masses, at the 3-0" confidence level, 
assuming a circular orbit in Figure [TU] We make use 
of the expre ssions for the TT V, TDV-V and TDV-TIP 
presneted in iKippingI ()2009al lbl) . We find that Kepler is 
clearly se nsitive to sub- Earth mass exomoons, as pre- 
dicted bv lKippind (|2009t ). 



Figure [TO] shows that for moons co-aligned to the 
planet's orbital plane {is — 90°), moons down to sub- 
Earth mass are excluded. The sensitivity drops off 
as inclination increases away from a co-aligned system 
but stabilizes for highly inclined moons (the kinks close 
is 0° and is ^ 180°) as a result of TDV-TIP effects 
dominating. 

6.5. Proposed Inclination Change 

iMislis fc Schmitd ()2009| ) claimed to have detected a 
linear decrease in the duration of TrES-2b due to the 
inclination angle varying at a rate of —0.195° over 
^300 cycles, or —0.00065° per cycle. 

Because the other ground-based measurements did not 
have their limb darkening coefficients fitted for, a fair 
comparison is not possible, in our view. Although the 
expression for the duration is independent of limb dark- 
ening parameters, we found that the duration was highly 
correlated to the limb darkening coefficients. However, 
we are able to use our 18 measurements of the inclina- 
tion to quantify the constraints on the rate of inclination 
change in this system. Comparing data taken from the 
same instrument which has a constant CCD response 
function and bandpass is justifiable without fitting for 
limb darkening each time, since the LD parameters will 
not vary transit-to-transit. 

Fitting a linear trend through our inclination data 
gives a rate of change of -f (0.0019 ± 0.0020)° per cy- 
cle, which is clearly not significant. We exclude an incli- 
nation change of —0.0041° per cycle to 3-(t confidence, 
which is larger than that claimed by [Mislis fc SchmittI 
(|200l . Therefore, using the current Kepler data alone 
is not sufficient to yet confirm/reject the proposed in- 
clination change in this system, largely due to the very 
small temporal baseline of just 18 cycles. We will return 
to this hypothesis in our study on the long term timing 
changes in ^7.2\ 

6.6. Other Changes 

6.6.1. Baseline 

The baseline fluxes are in excellent agreement with the 
global mean at all epochs, giving = 18.8 for 16 degrees 
of freedom. This is not a surprise since the baseline has 
been normalized twice during our corrective procedure 
(see !j2.2|) to ensure precisely this result. 

6.6.2. Transit depth 

The transit depths are extremely stable yielding = 
8.7 for 16 degrees of freedom, which is our most sta- 
ble statistic. The TSVs (transit depth variations) are 
shown in Figure llll where the low scatter, of standard 
deviation 59.3 ppm, is evident. We exclude variations of 
123 ppm to 3-(T confidence. Over the timescale of years, 
transiting planets on periods > 10 days may exhibit vari- 
ation s due to precession o f an oblate planet's rotation 
axis (jCarter fc WinnI l2010 f). However, we do not pos- 
sess sufficient baseline to look for such effects with the 
18 cycles of this study. 

7. LONG-TERM TIMING VARIATIONS 

7.1. Ephemeris Fitting 

In Table [5] of the appendix, we show all of the mea- 
surements of the transit times of TrES-2b used in this 
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Fig. 9. — TDVs of TrES-2b. Top left: Observed TDVs for TrE3-2b using the duration of the global fit. Bottom left: F-test periodogram 
of the observed TDVs. Top right: TDVs computed from control data (artificial lightcurves). Bottom right: F-test periodogram of the 
control TDVs. 
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Fig. 10. — Excluded exomoon masses for a companion to TrES-2b, as a function of the orbital distance of the moon around TrES-2b 
and the orbital inclination with respect to the orbital plane of TrES-2b. Contours are given in units of Earth masses, making steps of 
0.25 Mq. The Kepler data is able to easily probe down to sub-Earth mass exomoons. 



study, including both amateur and professional mea- 
surements. The inclusion of the previous transit times 
leads to much tighter constraints on the period and 
epoch. We repeated our fits without using the pre- 
vious transit times and found a local period of P = 
2.4706112tg °°™°^^days. Using all of the transit times 
yields P = 2.470618921°;°™°™^^ days, which is slightly 
longer than that found using the Kepler data only (note 



the much higher precision of using all of the transit 
times). This discrepancy is visible in Figure [8] where 
a drift in the TTVs is apparent and an excess of low- 
frequency power exists in the periodogram. Whilst both 
values are consistent with the iHolman et al.l ()2007l ) val- 
ues of P = 2.470621 ± 0.000017days, the reason for this 
discrepancy warrants further investigation. 

The previous transit times have several differences to 
the Kepler times; they are mostly from amateur as- 
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Fig. 11.— Transit depth variations (T5V) of TrES-2b. 
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Fig. 12. — Marginalized posterior distribution of the rate of 
change of the orbital period for TrES-2b (C). We find no evidence 
for a long-term change in the planet's orbital period. 



tronomers and they have a longer temporal baseline by 
a factor of 31.5. Another difference is that the DAWG 
do not recommend basing scientific conclusions on the 
Kepler times to an absolute accuracy of less than 6.5 s, 
until such a time as this level of accuracy can be verified 
(see ). 6.5 s is certainly sufficient to explain the observed 
low-frequency power observed in Figure [H We therefore 
consider three possible hypotheses to explain the discrep- 
ancy: 

1. The amateur transit times are unreliable and bias 
our results 

2. There exists a long-term deviation away from a lin- 
ear ephemeris 

3. Systematic error in the Kepler times 

7.1.1. Hypothesis 1 - The ETD measurements are 
unreliable 

The list of transit times used in this study consists of 62 
amateur measurements from the ETD and 22 from peer- 
reviewed publications. The professional times should be 
considered reliable by virtue of the peer review process, 
however the amateur times may or may not be reliable. 
For several epochs, there are simultaneous measurements 
from both camps and one may use these to evaluate the 
reliability of the amateur measurements. 

Epoch 142 has two measurements from each camp (a 
total of four transit times). The weighted average o f the 
professional measurements from lRabus et all (|2009l) and 
iRaetz ct al. (2009) average to BJDtdb 2454308.46163. 
Each timing measurement d eviates frorn this point 
by 0.76-O- and 2. 18-a for the iRabus et al] (|2009[ ) and 
iRaetz et al.l ()2009() measurements respectively. The ETD 
amateur measurements deviate from this same point by 
0.13-cr and 0.76-0-, indicative of a highly consistent result. 

Epoch 278 has one fro m each camp , the professional 
measurement being from IRaetz et al.l (|2009f ) . The am- 
ateur measurement deviates from this point by 0.93-(t, 
even when negating the error on the professional mea- 
surement (0.74-(T when including both errors). 

Epoch 316 has one from eac h camp , the professional 
measurement being from Raetz et al.l ()2009f ) . The am- 
ateur measurement deviates from this point by 1.89-f7, 



even when negating the error on the professional mea- 
surement (Q.78-cr when including both errors). 

Epoch 395 has one from each camp, the professional 
measurement being from Mislis et alj (j2010[ ). The am- 
ateur measurement deviates from this point by 0.083-(T, 
even when negating the error on the professional mea- 
surement (Q.076-(T when including both errors). 

Finally, epoch 414 is contemporaneous with one our 
Kepler transits (Kepler epoch 10), t hree ETD measure - 
ments and one professional time from lMislis et al.l (|2010( ) . 
Neglecti ng the much snialle r error on our Kepler transit 
time, the lMislis et al.l ()201Clt ) time deviates away by 1.88- 
a. The ETD measurements deviate by 3.18-(T, 0.43-cr 
and 0.22-cr. 

In conclusion, the contemporaneous sample of eight 
amateur transit times indicates that the amateur mea- 
surements are highly consistent with both the profes- 
sional data and the Kepler times. Seven out of eight 
measurements were less than 1-cr away from the profes- 
sional determination, with one outlier at Z-a. This out- 
lier point would be disregarded automatically by our fit- 
ting algorithm anyway as a result of using median statis- 
tics (see i j3.1.3l) . Although we cannot perform this test 
on every signle amateur transit, it seems reasonable that 
the selected sample is a fair representation of the ETD 
database. 

7.1.2. Hypothesis 2 - Long-term non-linear ephemeris 

Having shown that hypothesis 1 is not supported by 
the current body of evidence, we move on to investigate 
our second hypothesis. We tried fitting all of the data 
again using a quadratic ephemeris throu gh the transit 
times, using the same approach as that of iHolman et al.] 
(jlOlQ). The model is tjm = tq -\- nP -\- n^C, where n is 
the epoch number and C is the curl. For a period which 
decreases over time, one expects C < 0. 

Using both the and statistics, we found that there 
is no preference for a quadratic trend in the transit times. 
In Figure I12[ we show the marginalized posterior distri- 
bution of C, which is symmetric about zero. The data 
excludes a change in the orbital period of 0.11 seconds 
per year, to 3-cr confidence. 

7.1.3. Hypothesis 3 - Systematic error in the Kepler 

times 
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The only contemporaneous transit measurement be- 
tween Kepler and the professional measurements is in- 
sufficient to test this hypothesis . This is b e cause the 
transit is question, measured by iMishs et al.l (|2010[ ) for 
epoch 414, has a precision of 518 s and is therefore not 
useful for testing Kepler^ a timing accuracy at the < 6.5 s 
level. 

In conclusion, the current body of evidence is insuffi- 
cient to determine the cause of the discrepant periods. 
However, further transits from Kepler will resolve this 
issue in the future. 

7.2. Duration Change 
7.2.1. Comparing different bandpasses 

As discussed in g631 iMislis fc SchmittI ()2009l ) claimed 
to have detected a linear decrease in the duration of 
TrES-2b due to the inclination angle varying at a rate 
of 0.195° over ~300 cycles, or 0.00065° per cycle. 

In order to look for evidence of long duration change, 
it is n ecessary to use data taken before the Kepler Mis- 
sion. iHolman et al.l ()2007() obtained three high-quality 
transits observations using the FLWO 1.2 m telescope, 
in anticipation of this requirement, with a mean cycle 
value of 20.6. In contrast, the QO and Ql data have a 
mean cycle value of 412.5, giving a baseline to the FLWO 
data of 391.83 cycles. 

However, as also discussed in §6.5) we cannot com- 
pare data from different bandpasses unless we fit for limb 
darkening coefficients, especially since transit parameters 
are known to be acutely correlated to lim b darkening for 
near-grazing transits (see H3.2.2l and ii4.ll) . iHolman et al.l 
([2007) did not fit for limb darkening coefficients and so 
we here present a re-analysis of those three transits. 

7.2.2. Re-analysis of \Holman et al\ \200^ ) photometry 

We choose to fit for linear limb darkening due to the 
lower signal-to-noise from ground-based data. U sing the 
corrected photometry from lHolman et aD ()2007| ). we per- 
form these fits in the same manner used in this paper. 
We float esinw, ecosw, K , 7 and P around their best-fit 
values from the global eccentric run (see Tabled]) to allow 
their errors to propagate into the MCMC. The results are 
reported in Table 01 

7.2.3. Choosing a statistic 

In ^6.11 we compared durations within the same band- 
pass and thus fixing limb darkening was justified. Since 
Ti is highly correlated to the limb darkening coefflcients, 
we find TiA offers the highest signal-to-noise for TrES-2b 
when limb darkening is fitted and will be adopted here. 

7.2.4. Limits on duration change 
The IHolman et all (j2007D global fit finds Ti^^ = 
6559^]^Q3S whereas the Kepler global fit finds Ti^4 = 
6439^28 8, giving ATi,4 = (-120 ± 106) s over 
391.83 cycles, which we do not consider to be signifi- 
cant. The data exclude a decrease in the transit duration 
|ATi^4| > 4 38s, or 165s per year , to 3-cr confidence. In 
contrast, Mis lis fc SchmittI (j2009( ) claim to have detected 
a duration decrease of ~ 3.16 minutes (189.6 s) over ~ 
300 cycles (or 91 s per year). Whilst this is not supported 
by our analysis, it is also not excluded. IScuderi et al.l 



TABLE 4 

Resu lts from re-analysis of TrES-2b using the \Holman et oil 
1200%) photometry. Quoted values are medians of MCMC trials 
with errors given by 1-cr quantiles. * = fixed parameter; f = 
parameter was floated but not fitted. 
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(201(f) have challenged the 'Mis lis fc Schm itt' (2009^ re- 
sult rece ntlty using ground-ba sed transit observations 
and the iGilliland et al.l (|2010t ) result. Future Kepler 
transits will resolve this issue definitively. 

8. SUMMARY OF RESULTS 

Due to the large number of results presented in this 
paper, we summarize the key findings below: 

■ The Kepler SC data exhibit unprecedented preci- 
sion with r.m.s. noise 237.2 ppm per 58.8 s. 

■ Fitting for limb darkening coefficients leads to 
much larger uncertainties in the system parame- 
ters of TrES-2b, due to the near-grazing nature of 
the orbit (e.g. a factor of 17.5 larger for the transit 
depth) 

■ We present a self-consistent, refined set of tran- 
sit, radial velocity and physical parameters for the 
TrES-2b system, which are in close agreement with 
previous values. 

■ We do not detect an occultation of TrES-2b, con- 
straining the depth to be < 72.9 ppm to 3-a con- 
fidence, indicating that this object has the low- 
est measured geometric albedo for an exoplanet, 
of Ag < 0.146. 

■ We detect no short or long term transit timing vari- 
ations (TTV) in the TrES-2b system and exclude 
short-term signals of r.m.s. > 7.11s and a long- 
term variation of 0.11s per year in the orbital pe- 
riod, to 3-CT confidence. 

■ We detect no short or long term transit dura- 
tion variations (TDV) in the TrES-2b system and 
exclude short-term relative duration change of > 
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0.77% and long-term change of > 2.5% per year, 
to 3-CT confidence. 

■ We exclude the presence of exomoons down to sub- 
Earth masses for TrES-2b. 

■ The iMislis fc SchmittI ()2009D hypothesis of long- 
term duration change is neither supported nor re- 
futed by our analysis. 

■ We find the iRabus et all ()2009f ) hypothesis of a 
0.2 cycle TTV is not supported by the Kepler pho- 
tometry, to a high confidence level. 

■ We find no evi dence for oth e r dips in the lightcurve 
as reported bv lRaetz et al.l ()2009t l. 
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TABLE 5 

All measured transit times of Tr ES-2b taken from the lite ratur e, the ETD (Exopla net Trasnit Database) and this work, at the time of 

writing. * = value presented in \Mislis & Schmiti S20091} and \Mislis et ali 120 1(X) for this transit do not agree with each other and 
therefore these measurements are not included in our long-term TTV analysis, f = value comes from summing more than one transit 
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